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Abstract 

A theory based on Manning roughness equation, Philip equation and water balance equation was developed which only 
employed the advance distance in the calculation of the infiltration parameters and irrigation coefficients in both the border 
irrigation and the surge irrigation. The improved procedure was validated with both the border irrigation and surge 
irrigation experiments. The main results are shown as follows. Infiltration parameters of the Philip equation could be 
calculated accurately only using water advance distance in the irrigation process comparing to the experimental data. With 
the calculated parameters and the water balance equation, the irrigation coefficients were also estimated. The water 
advance velocity should be measured at about 0.5 m to 1.0 m far from the water advance in the experimental corn fields. 
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Introduction 

Border irrigation and surge flow irrigation method have been 
applied in arid and semi-arid regions to conserve water and to 
increase water productivity since 1987 in China [1]. Both of these 
irrigation methods have the advantages of achieving more uniform 
irrigation, saving water resources by speeding up irrigation water 
advance and reducing deep percolation [2-3], which was wide 
spread in China. Infiltration properties of the soil and irrigation 
coefficients are among the most important factors in the design 
and management of irrigation systems, but it is often difficult to be 
obtained because of the spatial and temporal variability of soil 
properties [4-5] . Due to the importance of irrigation parameters 
in surface irrigation system design and evaluation, numerous 
investigations are being conducted to measure or predict the 
surface irrigation events such as advance, storage, recession and 
subsurface [6-10]. Even though direct measurements are techni- 
cally preferable, the parameters obtained may not be applicable at 
field scale [1 1]. As a result, indirect evaluation of field parameters 
by numerical inversion of the governing equations has become an 
effective alternative to directive methods for the estimation of 
surface irrigation parameters. Three are several methods to 
estimate the surface irrigation parameters, such as one-point 
method, two-point method, as well as advance-repercussion 
method. The one-point method and two-point method could only 
estimate the surface irrigation parameters in a flat, homogenous 
field [12-13]; moreover, the advance-repercussion method could 
obtain the surface irrigation parameter only with the mean water 
depth [14] or the water storage in the field after the irrigation [15]. 
The review of literature suggests that even though various methods 



are available for the estimation of surface irrigation parameters, 
the needed information, such as irrigation advance, recession, and 
irrigation water depth data for the parameters, are very hard to be 
obtained in the field with an accuracy value except of the advance 
time and water advance distance. 

Thus, the objective of this study was to develop a simple method 
based on Manning roughness equation, Philip equation and the 
volume balance equation to predicate the border irrigation and 
surge flow irrigation parameters and infiltration parameters only 
with the water advance trajectory which is easy to be measured 
accurately in a surface irrigation experiment. 

Model Development 

Border irrigation 

The unit width flux and water flow velocity for irrigation could 
be described by the manning equation as: 

n 

(2) 

n 

where, v is the water advance velocity, L/T; q is the unit width 
flux, V/T; h is the depth of flow, H; J is the slope degree; n is the 
manning roughness coefficient. 
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In order to analyze the irrigation advance process, we assume 
that the depth of irrigation water among the distance be described 
with an exponential function as follows: 

h=h 0 (i--r (3) 

Where h 0 is the depth of surface water at the inlet, H; x is the 
distance to the inlet, L; Xf is the surface water advance distance, L; 
m is a shape coefficient of the surface water profile. 

Furthermore, in order to determine the water advance velocity, 
we assume that the water advance velocity at x r a may represent 
the velocity of water advance velocity, and the corresponding 
water depth will be: 



equations is the Philip equation: 

I = St 05 (11) 

Where I is the cumulative infiltration, L; S is the sortivity, 
L/T 1/2 . 

The water starts infiltrating into the subsurface soil at a point 
along the border strip only after the irrigation front reaches that 
point. Denoting t x as the net infiltration time at V point, 

And then, 

tx = t Xf -L (12) 



h x = ho(-T 
x f 



Correspondingly, the water advance velocity at x r a is: 

2m/3 



(4) 



x f 



\ x f) 



(5) 



There is a functional relationship of the advance distance with 
the advance velocity: 



dxf 

-di =V *f 



(6) 



Connecting Eq. (5) with Eq. (6): 



Where, t x f, x t - are the time that the irrigation front reaches to x f 
point and x point, respectively. 

According to eq (1 1), the infiltration at the x f point could be 
described as: 



T _ 0,0.5 _ g , D _ 'D,' 



DsO.5 



(13) 



As a result, the total infiltration amount (F), from the inlet to the 
distance x f is: 



x f 



Idx 
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Integrating Eq. (7): 



(l+2w/3)-7°-V /3 (0) 2m/3 



l+2m/3 

K f 



(7) 



(8) 



Eq (18) could be simplified as follows: 



(l+2D)S i +D/ 2 
2(l+D)D^v° 0 5 a( D -V/ 2X f 



(15) 



When water reaches the point, x f , the mount of surface water 
storage is: 



Let D = l+2m/3, and then Eq. (8) becomes: 
n 



(9) 



Jo Xf l+m 



(16) 



Therefore, the water balance equation could be written as: 



and Eq. (8) could be written as: 



t = 



Dv 0 a L 



(10) 



Eq. (10) describe the relationship between the surface water 
advance distance and the time. 

Soil infiltration characteristics are usually expressed in a time- 
dependant infiltration equation. One of the most common 



qot -- 



(1+2Z>)S 



x f +D ' 2 +- 



h 0 x f 



2(l+D)D°- 5 v 0 0 5 a( D -^ 2 f l+3(D-l)/2 



(17) 



Surge irrigation 

Similar to the theory improvement in Border irrigation, the 
theory was improved in the same way. During the surge irrigation 
process, when each surge was conducted, the irrigation time was 
measured as t(i), 
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Table 1. Soil characteristics of the experiment field. 



Depth (cm) 


Density (g/cm 3 ) 


clay (%) <0.002 mm 


silt (%) 0.002-0.02 mm 


sand (%) 0.02~2 mm 


0-10 


1.53 


5.3 


22.7 


72 


10-20 


1.51 


6.6 


26.9 


66.5 


20-30 


1.48 


5.9 


24.8 


69.3 


30-40 


1.49 


5.9 


24.5 


69.6 


40-50 


1.45 


5.6 


23.3 


71.1 


50-60 


1.59 


2.4 


8.8 


88.8 


Average 


1.51 


5.3 


21.8 


72.9 


doi:1 0.1 371 /journal.pone.01 01 274.t001 



The unit width flux and water flow velocity for each irrigation 
surge could be described by the manning equation as: 



n(i) 



(18) 



(19) 



where, v(i) is the water advance velocity for each surge, L/T; q(i) is 
the unit width flux for each surge, V/T; h(i) is the depth of flow in 
each surge, H; J is the slope degree; n(i) is the manning roughness 
coefficient for each surge. 

In order to analyze the irrigation advance process, we assume 
that the depth of irrigation water among the distance be described 
with an exponential function as follows: 



h(i) = h 0 (i)(\- 



*(0 yn(i) 

x f (i) 



(20) 
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Figure 1. Comparison of the measured and fitted water advance process. 
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Table 2. Fitted parameters from the measured water advance process. 







replications 


V„xa A (1-D) 


D/m/s 1 ' 2 


R 2 


1 


2.006 


1.619 


0.998 


2 


2.239 


1.511 


0.986 


3 


1.636 


1.634 


0.998 



doi:1 0.1 371 /journal.pone.01 01 274.t002 
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Figure 2. comparing the measured infiltration curve of the border irrigation with the simulated infiltration curve. 

doi:1 0.1 371 /journal.pone.01 01 274.g002 



Where h 0 (i) is the depth of surface water at the inlet in each 
surge, H; x(i) is the distance to the inlet in each surge, L; Xf(i) is the 
surface water advance distance in each surge, L; m(i) is a shape 
coefficient of the surface water profile in each surge. 

Furthermore, in order to determine the water advance velocity, 
we assume that the water advance velocity at Xf^-a(i) may 
represent the velocity of water advance, and the corresponding 
water depth will be: 



«(0 xm© 



Xf(l) 



(21) 



and Eq.(8) could be written as: 



t(t) = ^ — r x(i)f 



(27) 



Eq.(27) describe the relationship between the surface water 
advance distance and the time. 

Soil infiltration characteristics are usually expressed in a time- 
dependant infiltration equation. One of the most common 
equations is the Philip equation shown as follows, 



Correspondingly, the water advance velocity V x (i) at Xf(i)-a(i) is 



■n 

2m(/)/3 



7(0 = 5(0/(0° 



(28) 



1 / a(i) \ an W i Where I(i) is the cumulative infiltration in each surge, L; S is the 

v Y (0 = ^j jl/2h °(') 2P (t^) ( 22 ) sortivity in each surge, L/T 1/2 . 

The water starts infiltrating into the subsurface soil at a point 
along the border strip only after the irrigation front reaches that 
There is a functional relationship of the advance distance with P oint Denoting t x (i) as the net infiltration time at 'x' point, 



the advance velocity: 



And then, 



dxf(J) 

-dm =v *f {1) 



(23) 



t x (i) = t Xf (i)-t x (i) 



(29) 



Connecting Eq. (22) with Eq. (23): 



dt(i) nil) Wo \x f (i)) 



Integrating Eq.(24) 



i 



(24) 



( 1 + 2m©/3) / a5 /*o(0 2/3 (a(0) 2 '" ( ' )/3 f 
n(i) 



,l+2m(0/3 



(25) 



Where, txf, x t > are the time that the irrigation front reaches to Xf 
point and x point in each surge, respectively. 

According to eq (29), the infiltration at the x f point in each surge 
could be described as: 



/(0 = ^(0?(0°' 5 = 
S(i) 



, 0 Vn 0 ^n< fl «-i)/2 



'BnO.5 



£(0 U3 v(OoXO' 



(30) 



As a result, the total infiltration amount (F), from the inlet to the 
distance x f is: 



Let D(i) = l+2m(i)/3, and then Eq.(8) becomes: 



v 0 (i)=—h 0 (if 3 J 1 / 2 
n(i) 



(26) 



Table 3. Fitted parameters with the improved theory. 





replications 


aim 


A m/s A0 5 


h 0 /cm 


n 


1 


1.003 


0.048 


3.9 


0.036 


2 


1.200 


0.073 


4.1 


0,038 


3 


0.818 


0.049 


3.7 


0.035 


average 


1.007 


0.057 


3.9 


0.036 



doi:1 0.1 371 /journal.pone.01 01 274.t003 
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Figure 3. Water advance process in three cycles of Surge irrigation. 

doi:1 0.1 371 /journal.pone.01 01 274.g003 
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Eq (31) could be simplified as follows: 
(1 +2D(0),S(0 



m- 



2{\+D(i))D(i)^v 0 {i) u 'a{i) 



Xf{i) 



1 +23(0/2 



(32) 



•*/(0 



^o(0d -*«r ( 'V*=M^ (33) 



Therefore, the water balance equation could be written as: 
q 0 (i)t(i) = 

(1+2Z?(Q)S(Q m l+D»/2 

2(l+X>(0)X>(O a5 Vo(O a5 «(0 (i ' (, ' )_1)/2 (34) 



+ 



h 0 {i)xf(i) 



l + 3(/)(0-l)/2 



When water reaches the point, Xf(i), the mount of surface water 
storage in each surge could be described as, 



Materials and Methods 

The field experiments were conducted in Bayinguoleng 
Mongolia Autonomous prefecture Bazhou irrigation station which 
was one of the field stations of our Scientific team (Prof. Wang 



Table 4. Mean fitted 


parameters of the water advance process 


in all the three cycles. 








D m/sl/2 


V 0 xa A (1-D) 


R 2 


First cycle 


1.6771 


1 .4791 


0.9931 


Second cycle 


1 .4926 


1.5053 


0.9694 


Third cycle 


1 .6984 


0.4236 


0.9876 



doi:1 0.1 371 /journal.pone.01 01 274.t004 
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Figure 4. Estimated infiltration curve of the border irrigation. 
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Quanjiu was response for experimental fields in this manuscript 
and can be contacted in the future). And I declare that my 
manuscript does not matter the any ethic issue. 

Bayinguoleng Mongolia Autonomous prefecture Bazhou irriga- 
tion station was classified as a warm-temperature arid zone with 
continental climate. The field experiments were carried out on 3 
bays (4*200 m) in a corn field. Corn is an important crop here 
which requires an average of 6 irrigations; however, the irrigation 
time may be increased to 7-9 irrigations in drier climate. 

The measured values of the physical properties of the soil are 
given in Table 1 , including the soil bulk density and percentage of 
sand, silt and clay, at the depths of 0-10 cm, 10-20 cm, 20- 
30 cm, 30-40 cm, 40-50 cm, 50-60 cm. The soil at the 
experimental site is classified as sandy soil. 

In both of the border irrigation and surge irrigation experi- 
ments, irrigation water were pumped onto the bay through a weir 
plate which could measure the water flow rate, and then was shut 
off when the advance front reached a distance of 200 m. The flow 



rate was about 3.94 L/m-s. Border irrigation advance time was 
recorded that when the advance front arrived at a known distance. 
The soil water contents of the soil profile before and after the 
border irrigation were measured and then be used to verify our 
improved theories. In our surge irrigation, there were three cycles 
of on-time and off-time periods. Water advances to the end of the 
field only during the last surge. During each surge, water advances 
and infiltrates in the same way as that in a continuous irrigation. In 
each cycle, the water contents in the soil profiles before and after 
the surge irrigation were also measured, as well as the infiltration 
rates measured with a two-double ring method in order to certify 
the proposed theory. Both the border irrigations and surge 
irrigations were conducted 3 times. 

Results and Discussions 

Border irrigation 

Trials of the models have been carried out using three 
experimental data sets for discussions. In order to understand 



Table 5. Mean fitted parameters in three irrigation cycles. 






First cycle 


Second cycle 


Third cycle 


A/m/s A °- 5 


0.081 


0.029 


0.0018 


H 0 /cm 


3.7 


2.0 


2.0 


N 


0.020 


0.007 


0.007 


v 0 /m/s 


0.35 


0.639 


0.639 


a/m 


1.23 


0.587 


0.586 


doi:1 0.1 371 /journal.pone.01 01 274.W05 
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the water advance phase under the border irrigation, results from 
experimental site were shown in Fig. 1 . From Fig. 1 it can be easily 
found that the water advanced distance increased with the time 
increased. In order to evaluate the theory, the water advancing 
distance changing with the time was fitted well by a power 
function and the parameters were shown in Table. 2. 

Based on the measured water advanced curves, the parameters 
of VoXa A (l-D) and D were easily obtained, thus the parameters of 
S, a, h 0 could be obtained with eq(10) and eq(17) according to the 
volume balance equation. Once the parameter of h 0 is known, the 
roughness coefficient could be calculated with equation (1) or 
equation (2) shown in Table 2. 

From Table 3, we can found that the average value of 
parameter a may represent the water advance velocity before 
the water advance 1.007 m that may represent the true water 
advance velocity. The irrigation water depth was about 3.9 cm 
which was also credible in the irrigation event. Then with the 
parameters above, the cumulative infiltration changing with the 
time could be estimated as follows shown in Fig.2. From Fig.2, it 
could be found that the new developed theory could estimate the 
infiltration process very well. 

Surge irrigation 

Similar to Border irrigation process, in order to certify our 
improved theory about the surge irrigation, the surge irrigation 
experiments were also conducted in the Corn field. The 
experimental field was 2.5 m in width, 200 m in length. Each 
experiment includes three cycles of on-time and off-time periods. 
The time for the first surge irrigation starts at 0 min, the second 
surge irrigation starts at 120 min, and the last surge irrigation 
starts at 240 min. Water advances to the end of the field only 
during the last surge. During each surge, water advances and 
infiltrates in the same way as that in a continuous irrigation. 

In order to describe the surge irrigation process clearly, the 
advance distance in a certain time was shown in Fig. 3 and the 
process was fitted by a power function. Then the needed 
parameters could be obtained from the fitted results shown in 
Table 3. 

According to the equation(19) to (34), the parameters could be 
calculated as follows in Table 4. 

Table 5 shows the mean fitted parameters in three irrigation 
cycles. The infiltration processes of all the three surge irrigations 
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